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ABSTRACT 

We present a new IRAC, Spitzer IRAC images of the HH 34 outflow. These are the first images that 
detect both the knots along the southern jet and the northern counterjet (the counterjet knots were 
only detected previously in a long slit spectrum). This result removes the problem of the apparent 
coexistence of a large scale symmetry (at distances of up to ~ 1 pc) and a complete lack of symmetry 
close to the source (at distances of ~ 10 17 cm) for this outflow. We present a quantitative evaluation 
of the newly found symmetry between the HH 34 jet and counterjet, and show that the observed 
degree of symmetry implies that the jet production region has a characteristic size < 2.8 AU. This is 
the strongest constraint yet derived for the size of the region in which HH jets are produced. 
Subject headings: circumstellar matter — stars: formation — ISM: jets and outflows — infrared: ISM 
— Herbig-Haro objects — ISM: individual objects (HH34) 



1. INTRODUCTION 

HH 34 is one of the HH objects in Herbig's catalogue 
(Herbig 1974). It jumped into the limelight with the pa- 
per of Reipurth et al. (1986), who showed that the HH 
object (HH 34S) had a bow shaped morphology, and a 
jet-like association of aligned knots (pointing towards the 
apex of the bow shock). Later observations showed the 
existence of a Northern counterpart to HH 34S (HH34N) 
and of a series of bipolar bow shock pairs at larger dis- 
tances from the outflow source (Bally & Devine 1994; 
Eisloffcl & Mundt 1997; Devine et al. 1997). There is 
a wealth of observations of this outflow, including im- 
ages, spectrophotometry, radial velocity and proper mo- 
tion measurements at optical (see, e. g., Heathcote & 
Reipurth 1992; Eisloffel & Mundt 1992; Morse et al. 
1992, 1993; Reipurth et al. 2002; Beck et al. 2007) 
and IR wavelemgths (Stapelfeldt et al. 1991; Stanke et 
al. 1998; Reipurth et al. 2000). These observations show 
that the HH 34 outflow (at a distance of « 417 pc, see 
Menten et al. 2007) has a plane of the sky velocity of 
« 150 km s _1 and propagates at an angle of « 30° from 
the plane of the sky (the southern lobe being directed 
towards the observer). 

The flow itself is driven by HH 34 IRS, a Class I proto- 
star surrounded by a relatively large (~ 1000 AU radius) 
circumstellar disk (Stapelfeldt & Scoville 1993; Anglada 
et al. 1995). The possible coupling of the disk and out- 
flow has motivated the search for signatures of internal 
rotation within the jet as consequence of the transfer of 
angular momentum from the rotating disk into the highly 
collimated jet (Coffey et al. 2010). 

In spite of the clear bipolar symmetry of the bow shock 
pairs (extending to ~ 1.5 pc from the source, see Devine 
et al. 1997), observations extending over ~ 2 decades 

1 Institute de Ciencias Nucleares, Universidad Nacional 
Autonoma de Mexico, Ap. 70-543, 04510 D.F., Mexico 

2 SPITZER Science Center, California Institute of Technol- 
ogy, C A 91125 USA 

3 Astronomisches Rechen-Institut Zentrum fur Astronomie der 
Universitat Heidelberg, Monchhofstr. 12-14 69120 Heidelberg, 
Germany 

4 Jet propulsion Laboratory, California Institute of Technol- 
ogy.MS 183-900, 4800 Oak Grove Drive, Pasadena, CA 91109, USA 



did not detect a nothern counterpart for the chain of 
aligned knots extending ~ 30" southwards from the HH 
34 source. This situation changed with the paper of 
Garcia Lopez et al. (2008), who obtained IR (1.6 and 2.1 
(im) long-slit spectra in which the emission of the north- 
ern counter-jet was finally detected. These authors note 
that the emission of the counterjet has intensity peaks 
at positions (i. e., distances from the source) which ap- 
proximately coincide with the knots along the southern 
jet. 

In this paper, we present new Spitzer IRAC images of 
HH 34. These images show the southern jet and north- 
ern counterjets with comparable intensities, and with a 
surprising degree of symmetry. The observations are de- 
scribed in section 2. In section 3 we present an image of 
the central region of the HH 34 outflow, quantitatively 
evaluate the degree of symmetry between the jet and the 
counterjet, and discuss the implications of the results for 
the ejection mechanism that has produced the outflow. 
The results are summarized in section 4. 

2. OBSERVATIONS 

The observations of HH 34 are part of our original 
Spitzer Space Telescope (Werner et al. 2004) General 
Observer (GO) program 3315 (PI Noriega-Crespo) ob- 
tained with both the infrared camera IRAC (Fazio et al. 
2004) and the infrared photometer MIPS (Rieke et al. 
2004) in March 28, 2005. The data have been recovered 
from the Spitzer Legacy Archive and the quality of the 
final images (Post Basic Calibrated Data or Post-BCD; 
S18.7 products) is outstanding, so that no further pro- 
cessing was required. In this study we present the IRAC 
observations obtained in the four channels (1, 2, 3, 4) = 
(3.6, 4.5, 5.8 & 8.0 ^m) covering a FOV of - 30 x 30' (the 
result of a 6 x 6 array map with a 260" stepsize) and with 
a total integration time per pixel of 30 sec. The final 
images are sampled with 0.6" per pixel, nearly one third 
of standard ~ 2" IRAC angular resolution. 

Figure 1 shows a three color image of HH 34 using 
channels 1, 2 and 3. Like with other protostellar out- 
flows observed with IRAC (see e.g. Noriega-Crespo et al. 
2004; Looney, Tobin & Kwon 2007; Tobin et al. 2007; 
Ybarra & Lada 2009) Channel 2 recovers the strongest 
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jet emission, since its bandpass (4 to 5/Ltm) includes three 
relatively bright pure rotational H2 emission lines, 0-0 
S(9) 4.69, 0-0 S(10) 4.41 and 0-0 S(ll) 4.18/xm; the jet is 
well detected in Channel 3 as well, where another couple 
of H 2 lines are found (0-0 S(6) 6.11 and 0-0 S(7) 5.51>m) 
(see De Buizer & Vacca 2010). 
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Fig. 1— A section of the HH 34 IRAC map (~ 3'x 4') cen- 
tered on the HH 34 IRS source using the four channels 3.6 (blue) 
, 4.5+5.8 (green) and 8fim (red). North is up and East to the left. 
The image shows the counter-jet as well as a wider infrared bow- 
shock about 20" South of the optical HH 34N main bow shock. 
Indeed both HH 34S and 34N lie within the image, but are barely 
visible. Other extended emission structures (e. g., the North-South, 
jet-like structure South of HH 34 IRS) seen in the map might be 
associated with other outflows in the region. 



3. THE JET/COUNTERJET SYMMETRY 

In Figure 2, we show the 4.5+5.8/im image of the re- 
gion around the source of HH 34, rotated so that the out- 
flow axis is parallel to the ordinate. This image shows a 
surprising symmetry between the southern jet (detected 
in optical images) and the northern counterjet. We have 
defined a 7 pixel (4". 2) wide, recangular box aligned with 
the HH 34 axis, within which we subtract the background 
emission (defined as a linear interpolation between the 
pixels immediately outside of the long edges of the box) . 
The result of this background emission subtraction is 
shown on the right hand side of Figure 2, in which the 
knot structure along the HH 34 jet and counterjet is seen 
more clearly. 
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Fig. 2. — Central region of the 4.5+5.8 (im IRAC map (see Fig- 
ure 1), rotated so that the HH 34 axis is parallel to the ordinate. 
The images correspond to the addition of the two bands, and are 
displayed with the logarithmic colour scales given (in MJy/sterad) 
by the bars on the top of each plot. The left plot shows the map, 
and the right plot shows a map in which a background substrac- 
tion has been applied within a rectangular box that includes the 
jet and the counterjet (see section 3). The scale and orientation 
of the images are shown on the left hand plot. On the right hand 
plot, the red crosses show the positions of the outflow source and 
the jet /counterjet knots determined with the paraboloidal fits de- 
scribed in the text. The black circle shows the position of the source 
calculated as the average of the positions of the 7 jet/counterjet 
knot pairs (see the text). 



We have carried out paraboloidal fits to the emission 
peaks of the knots along the outflow axis, determining 
the positions of the peaks with an error of w 0.2 pix. 
= 0".12. Through this procedure, we obtain the posi- 
tions Xj and x c j of the 7 knots along the jet and the 
counterjet (respectively). The position of the source is 
not so well determined, because the image of the source is 
partially saturated, and has a complex point spread func- 
tion. Because of this, we have estimated the position of 
the source as the average of the coordinates of the seven 
knots along the jet and the counterjet. A paraboloidal 
fit to the emission of the region around the source actu- 
ally results in a similar position (with offsets of « 0.2 pix 
along and ks 1 pix across the jet axis with respect to the 
average position of the ensemble of knots) . 

In Figure 3, we plot the positions x C j of the consecutive 
knots along the counterjet as a function of the positions 
xj of the corresponding knots along the jet. This figure 
shows the remarkable symmetry (with respect to the po- 
sition of the source) of the knots along the jet and the 
counterjet. 

In Figure 3, we also plot the offsets Ax = \xj — x c j | (be- 
tween the positions of the corresponding jet/counterjet 
knot pairs) as a function of Xj . The first three knots have 
monotonically growing offsets, with Ax = 0" .11 — > 0".44 
(ranging from ~ 1 — > 4 times the measurement error, 
see above) . The four knots further away from the source 
generally have larger offsets, with a top value Ax = 1".57 
for the 6th knot out from the source. The resulting Ax 
vs. Xj dependence therefore has low offsets close to the 
source, and increasingly large (and more highly oscillat- 
ing) values for larger distances from the source. 

These results can be interpreted as follows. Let us first 
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Fig. 3. — The top graph shows the distance from the source 
x c j of the knots along the northern counterjet as a function of 
the positions Xj of the corresponding knots along the southern jet. 
The bottom plot shows the offsets Aa: = \xj — x c j\ in the positions 
of the corresponding jet/counterjet knots as a function of distance 
Xj from the source. The vertical size of the crosses in the Air vs. 
Xj plot approximately correspond to the error in the observational 
determinations of the Ax values. 



assume that the jet/counterjet knot pairs are ejected at 
the same time r, but with different velocities Vj and v c j = 
Vj — Av (for the jet and counterjet knots, respectively). If 
the knots are ballistic, at a time t they will have positions 
Xj = (t — r)vj and x c j = (t — r)(vj — Av). One then 
obtains that the jet/counterjet knot offset Ax v (due to 
asymmetries in the jet/counterjet ejection velocities) is 
given by : 

_ _Av 

* tJb u — j OC £j — — - <Xj j ■ ^ \~ J 

V 3 

In other words, an ejection velocity asymmetry in the 
simultaneous ejection of a jet/ counterjet knot pair results 
in a knot position asymmetry Ax v that increases linearly 
with distance from the source (as the knots travel away 
from the source) . 

Let us now assume that we have a pair of knots ejected 
with the same velocity Vj , but at times r and r + At (for 
the knot along the jet and the counterjet, respectively). 
For ballistic knots, at a time t they would then be at 
distances Xj = (t — r)vj and x C j = (t — t — Ar)vj from 
the source. The jet/counterjet knot offset Ax T (due to a 
time-defference At in the ejection time) is given by : 

Ax T = Xj — x c j = Arvj . (2) 

In other words, the position asymmetry Ax T does not 
change as the knots travel away from the source along 
the jet and counterjet. 

For the case of jet/counterjet knot pairs ejected with a 
time difference At and a velocity difference Av, we there- 
fore predict an asymmetry between the jet and counterjet 



knot positions of the form 

Av 

Ax = \Ax v + Ax T \ = | — Xj+AtvjI . (3) 

v i 

This relation gives the behaviour of Aa; as a function of 
distance from the source Xj (along the jet) for a given 
pair of knots as they travel down the jet/counterjet axis. 

If we have a series of knot pairs (7 knot pairs in our 
HH 34 image), the values of Av/vj and Atvj will not 
necessarily be identical for all of the knot pairs. We can 
anyway attempt to fit a linear Ax vs. Xj dependence to 
the observations. 

At large enough distances from the source, equation §5§ 
is a straight line which goes through the origin (of the Aa; 
vs. Xj coordinates), with a slope of \Av/Vj\. The Arvj 
term (see equation [3]) results in a scatter of points above 
and below (depending on the sign of At this straight 
line. This scatter will also have a contribution from the 
distribution of the |Av/i>j| values over the knot pairs. 

Therefore, if we have a distribution of At (over the dif- 
ferent knot pairs) with a zero mean velocity (i. e., with 
similar positive and negative At values), we would ex- 
pect to recover from the data a straight line that still goes 
through the origin. The slope of this line would repre- 
sent an average of | Av/vj | over the ensemble of knot pairs 
(see equation [3]) . The value of the determined error for 
the Xj = interception point of the straight line would 
give us an upper limit to the root mean square average 
of Arvj. 

Fitting a least squares linear dependence to the ob- 
served offsets, we obtain 

Vj At = (-0.04 ±0.36)"; — = 0.038 ± 0.019 . (4) 

V 3 

We therefore obtain an estimate for the average time 
offset of basically (as expected), an upper limit for 
its root mean square value of 0".36/uj, and a velocity 
asymmetry with an average modulus of approximately 
(0.04 ±0.02)«j. 

If we consider a distance of 417 pc (to HH 34) and 
that the knots have a plane of the sky velocity Vj « 
150 km s , we then obtain 

AT<4.5yr; Av = (5.7 ± 2.8) kms" 1 . (5) 

The maximum value of At that we have determined in 
this way from the observed jet/counterjet asymmetry of 
HH 34 directly implies a maximum possible size for the 
region in which the ejection is produced. The gas at the 
foot of a stellar jet is observed to have temperatures of ~ 
10 3 K, implying a sound speed of ps 3 km s _1 . The Alven 
velocity in this region is expected to have similar values. 
In order to synchronize the jet/counterjet ejections to a 
time At, the region producing both outflows has to be 
able to communicate causally within this time interval. 
Therefore, waves propagating at a velocity v s (this could 
be the Alfven or the sound velocity, both ~ 3 km s _1 , 
see above) should be able to travel the characteristic size 
D of the jet formation region within a time At. This 
condition can be written as : 

D<v s At, soD< ( - Vs ) 2.8 AU, (6) 
\3kms L / 

where for the second inequality we have used the value 
of At derived from the observed HH 34 jet/counterjet 
asymmetry (see equation [5]) and the estimated sound 
speed of 3 km s _1 . 
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4. SUMMARY AND CONCLUSIONS 

We have obtained IRAC images of the region around 
the source of the HH 34 outflow (see Figure 1). These im- 
ages show the presence of a southern jet and a northern 
counterjet (previously only detected in the IR long-slit 
spectrum presented by Garcia Lopez et al. 2008). The 
seven knots along the jet and the counterjet have a re- 
markable symmetry with respect to the position of the 
outflow source (see Figure 2). 

We find that if we plot the offsets Ax — \xj — x C j\ 
(where Xj is the distance from the source of a knot along 
the jet and x C j the distance to the corresponding knot 
along the counterjet) as a function of the position Xj of 
the knots, we obtain a general trend of lager Ax values 
for increasing distances from the source (see Figure 3). 
This trend can be interpreted (on the basis of a simple, 
ballistic knot model) as the result of a time-offset At 
(between the ejection of the two knots in a given pair) 
and of a velocity difference Av (between the velocities of 
the knot pair). 

From a fit to the observed Ax vs. x relationship, we de- 
termine a velocity asymmetry of Av ~ (5.7±2.8) km s , 
and an upper limit for the time-offset of the binary ejec- 
tion of At = 4.5 yr. Interestingly, the velocity asymme- 
try estimated for the binary ejections is within a factor 
of ~ 2 of the estimates for the sound and Alfven speeds 
at the base of an outflow from a young star. Actually, 
the estimates of the sound and Alfven speeds lie within 
the error bars of the Av value determined from the HH 
34 jet/counterjet knot asymmetries. 

It is possible to use the upper limit of 4.5 yr to 
derive an estimate D w 2.8 AU for the size of the jet 
production region. For example, if one considers mod- 
els of magnetocentrifugally driven disk winds, this esti- 
mate would imply a characteristic radius for the Alfven 
surface (in which the outflow becomes super- Alfvenic) 
of D/2 » 1.4 AU. This is an interesting constraint on 
models of the production of jets from young stars (for a 
recent discussion of jet production models see the review 



of Ferreira 2008, and the papers in the book of Garcia & 
Ferreira 2010). 

Another possible interpretation is that the » 1,4 AU 
radius (derived from the observations) corresponds to the 
outer radius of the region of the accretion disk which 
gives rise to the jet. This interpretation might be more 
comfortable in terms of magnetocentrifugally driven disk 
winds, because an Alfven surface of ~ 1 AU radius might 
be too small. 

VLA observations of sources of HH outflows had pro- 
vided a constraint of a maximum size of rs 50 AU (0".l 
at the distance of Orion) on the size of the jet production 
region (see, e. g., the HH 1/2 observations of Rodriguez 
et al. 2000). The upper limit of « 2.8 AU that we 
have derived above for the size of the formation region of 
the HH 34 jet is therefore the most stringent condition 
yet. This size corresponds to ~ 7 mas, a resolution that 
will only be approached in direct observations of an HH 
source by future efforts with the ALMA interferometer. 
Observations with large baseline near-IR interferometers 
(such as Keck-I or VLT-I) could possibly also resolve such 
a region, provided that the emission is strong enough to 
be detected. 

To conclude, we should point out that the derived up- 
per limit of f» 2.8 AU for the jet formation region is 
derived from a simple, ballistic approximation for the 
motion of the knots along the HH 34 jet and counterjet. 
Of course, the observed asymmetry could partially arise 
from deviations from a purely ballistic motion. If this 
is the case, the size implied for the jet formation region 
would be even smaller. 
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